CKM matrix unitarity and a novel type of global fits 
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The aim of the paper is to propose one paradigm change of CKM global fits on experimental data 
from electroweak sector. The change refers to using in fits the exact unitarity constraints expressed 
in terms of four invariant parameters, such as moduli of the CKM matrix, and to take into account 
an important set of the available experimental data. In the paper we use data from nuclear beta 
decays, and from leptonic and semileptonic decays, in order to find the most probable numerical 
form of the CKM matrix, as well as the determination of decay constants, fp, and of various form 
factors /+ p (0), directly from experimental results. 

PACS numbers: 12.15.-y, 12.15 Hh, 12.15 Ff 



The consistency problem of experimental data with 
unitarity constraints was recently solved, and a proce- 
dure for recovering the CKM matrix elements from error 
affected data was provided in These unitarity con- 
straints say that the four independent parameters Sij and 
cos 5, with the standard notation, should take physical 
values, i.e. Sij G (0, 1) and cos<5 € (— 1, 1), when they are 
obtained from the equations: 
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S 12 C 23 + S 13 S 23 C 12 + 2si 2 Si 3 S 2 3Cl2C23 COS 5, 

c 12 c 23 + S 12 S 13 S 23 ~ 2si 2 Si 3 S 23 Ci 2 C 23 COS S, (1) 
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s 13 c 12 c 23 + s 12 s 23 _ 2si 2 Si 3 S 2 3Cl 2 C 23 COS 0, 

S 12 S 13 C 23 + c 12 s 23 + 2si 2 Si 3 S 23 Ci 2 C 23 COS S 



where Vij = \Uij\, and UijSse the entries of the CKM 
matrix. It was shown in |l[ that if the physical quan- 
tities could also depend upon CKM matrix moduli, the 
reconstruction of a unitary matrix from such data is es- 
sentially unique, and in the following the used indepen- 
dent parameters will be the moduli Vij. For example, if 
V us = a, V u b — b, V cc i = d, and V c b — c is one set of four 
independent moduli, the relations ([T]) give the following 
solution 
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Because there are 58 groups of four independent moduli, 
one gets 165 different expressions for cos S and all have to 
take (roughly) the same value when computed from Eqs. 
([1]). Thus the x 2_ f un ction have to contain two kinds of 
terms: the first has to impose the fulfillment of unitarity 
constraints, and the second should take into account the 
physical quantities measured in experiments. 

Concerning the experimental data we will consider 
data coming from super-allowed + — > + nuclear beta 
decays, and from leptonic and semileptonic decays. 



In the standard model the decay rate for purely lep- 
tonic decay is given by 
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where Gf is the Fermi constant, Mp and me are the 
masses of the decaying meson, and, respectively, of the 
final lepton, U qq i is the corresponding CKM matrix ele- 
ment, and fp is the decay constant. In general one has 
also to take into account the radiative corrections which 
lead to a minor modification of the above formula. 

The physical observable for semileptonic decays, that 
depends on \U qq i\ and f(q 2 ), is the differential decay rate 
which under assumption of massless leptons is written as 



dT(H -> Piv t ) 
dq 2 
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where q = pn — pp is the transferred momentum, and 



\(q 2 ) = (M 2 H + M 2 
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is the usual triangle function, and f(q 2 ) is the global 
form factor which is a combination of f+(q 2 ) and f-(q 2 ). 
The experimenters provide numerical values for products 
of the form \U qq i /+(0)|, and in this paper we will use 
these numerical values. Hence from such experiments 
one measures, up to known factors, products of the form 



\U l3 \ 2 xf 2 , and/or, \U l} \ 2 x \f(q 2 )\ 2 



(7) 



It is clear that from such measurements one cannot find 
two unknown quantities, let's say, |/(<? 2 )| and \Uij\, if we 
have no supplementary constraints. Our point of view 
is that the unitarity constraints, which depend only on 
\Uij \ moduli, see relations ([I])-©, provide the necessary 
tool for the separation of moduli, and f(q 2 ), or fp. 

Before defining our type of global fit we make one natu- 
ral assumption, which is: the numerical values for all the 
measured moduli, |£7y|, must be the same irrespective of 
the physical processes used to determine them. The other 
parameters, such as the decaying constants fp, form fac- 
tors /+(0), gA/gv: etc., which parametrize the data from 
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each given experiment, are considered free parameters to 
be found from fit. 

The hrst piece containing unitarity constraints enter- 
ing the x 2 -function has the form 
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and the second component, which takes into account the 
experimental data, is 



X2 = > A — — 
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where c?i are the theoretical functions one wants to be 
found from fit, di is the numerical matrix that describes 
the corresponding experimental data, while a is the ma- 
trix of errors associated to d{. In the following our % 2 - 
function will be 
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Concerning experimental data we use the following. 
Knowledge on \U u d\ comes mainly from three different 
sources: a) super allowed, + — > + , nuclear beta decays, 
see (3, H, and 0], b) neutron beta decay, n — > pe^v, 
see [5|-[13|], c) and pion beta decay 7r + — ► 7r°e + ^, jl4| . 

The used data for the determination of the decay con- 
stants f n , fx, /b, /d, and f D + are from the papers 
d, QJ, DS, El, and, respectively, from [l|-[2|. Nu- 

are from the papers 
come from [37j and 
\U cd f^)llS K ^)U cs \ is given in [3£ 
tapers [41j-[50| provide data on |.F(1)?7 C &| 



merical results on \f+*(0)U u 
[H, those upon \f^(0)U ub 



and the ratio 
and Eel Thei 



and [42j and [5l|-[53| provide values for \G(l)U c b\. 

The central values and uncertainties used in fit are 
those published in the above papers, and we combined 
the statistical and systematic uncertainties in quadrature 
when experimenters provided both of them. 

According to fsj, the super-allowed beta decays be- 
tween T = 1 analog + states, together with the con- 
served vector current (CVC) hypothesis, lead to the con- 
clusion that the ft values should be the same irrespective 
of the nucleus, i.e. 
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where K is the vector coupling constant for semi-leptonic 
weak intercations, / is the statistical rate function, and 
t is the partial half-life. Because the above relation is 
only approximately satisfied, one defines a "corrected" 
Tt value, which should be "constant" , as 



Tt = ft(l-5 R )(l-5 c ) 
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where So is the isospin-symmetry-breaking correction, 5r 
is the transition-dependent part of the radiative correc- 
tion, and A)^ is the transition-independent part. Numer- 
ical values for Tt are given in 0, [U, and [H 
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In our fit we 

2 , by supposing 
v 



use the above formula with \Gy\ 

that <?y(0) = 1, as CVC requires, and \U u d\ and are 
the free parameters to be obtained from fit. Similarly for 
the neutron beta decay data we make use of the formula 



|tU 2 (l + 3A 2 ) 



4908.7(1.9) s 
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see [54j [. where r„ is the neutron mean life, and the free 
parameters are \U u d\ and A = gA/gv- 

If one or more of the above parameters could be mea- 
sured in other experiments, this approach allows us to 
take the results of these measurements into account. 
That is the case of the ratio gA/gv which enters in the 
measured asymmetry parameter Aq, see papers (55j-[58|. 
Their effect was a lowering of A to the value given in the 
Table, while by using only results from neutron beta de- 
cay data the value, A = -1.27092 ± 0.00394, is obtained. 

Values and corresponding uncertainties obtained from 
the fit are given in Table. 

The surprising result of our fit was that A^ is not 
transition-independent as it is usually assumed, see 
Refs. 01 -01 1 an d [59j]-[60]. F° r example, if one uses the 
data on nuclear beta decays from [3| , the A^ variation is 
from 2.193% for 22 Mg, to 2.579% for 54 Co nucleus. For 
this case the corresponding mean value and uncertainty 
are given in Table. 

If one makes use of Savard tt al. data, one gets 



A£ = (2.294 ±0.131) 



(14) 



and the values spreading is between 2.027% for 74 Rb, and 
2.429% for 34 C1. 

Our approach allows the use of all the seventeen values 
from 0, and one gets A^ = (2.362 ± 2.133)%, where 
the value of A^ for 42 Ti provided by the fit is negative 
A^ = —4.673%! However the result makes sense since 
the corresponding Tt is 3300 ± 1100 which is far away 
from the mean value given in which is around 3072. 
Hence the fit suggests us to throw out this value. By 
excluding also the 18 Ne and 30 S data, that lead to greater 
values than the mean by a factor of 3, and, respectively 
of 2, one obtains 



A R 



(2.364 ±0.182)% 



(15) 



Looking at the three Ap, values, that from the Table, 
and those obtained by using data from [3], Eq. (Tl4|) . and, 
respectively, from Eq. (TT5")) . one observes that they 
are compatible within the errors, and the better data are 
those coming from Ref. [3j. In all these three cases the 
errors provided by the fit are bigger than the theoretical 
estimates. We remind that the theoretical estimates are 



(old) 



A^(new) 



(2.40 ±0.08)%, and (16) 
(2.361 ±0.038)% (17) 
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given respectively, in Ref. [6(|, and 



Parameters 


Central Values and Errors 


Vud 


0.974022 ± 3.9 x 10~ 6 


V us 


0.226424 ± 3.9 x 10~ 6 


V ub 


(3.57604 ± 0.00002) x 10~ 3 


v cd 


0.226261 ± 3.9 x 10" 6 


v cs 


0.973324 ± 4.1 x 10~ 6 


v cb 


(38.0239 ±0.0002) x 10~ 3 


v td 


(9.28657 ±0.000035) x 10~ 3 


V ts 


(37.0454 ±0.0002) x lO" 3 


v tb 


0.999270 ±2.2 x 10~ 7 




(2.399 ±0.108)% 


9A/gv 


- 1.26924 ±0.00510 




(3. 104 ± 0.096)% 


u 


130.784 ± 1.323 


fx 


154.535 ± 1.990 


fx/ fir 


1.1816 ± 0.0272 


f B 


281.97 ±0.39 


f D+ 


220.1 ±0.8 




268.42 ± 11.22 


/**(0) 


956.8 ±9.1 


/f*(0) 


243.2 ±20.8 


f D, {Q)/f DK {0) 


0.833 ±0.006 


m 


941.9 ±78.7 


G(l) 


948.1 ± 149.8 



As one conclusion one can say that our approach does 
not confirm the (approximate) constancy of A^, and at 
the same time it allows a fine structure analysis of all 
nuclear beta decays, or, more precisely, of the present 
procedure for getting a "constant" Tt. The solving of the 
constancy problem of A)^ could require new ideas. One 
suggestion could be the use of the present approach, but 
now with a few more steps. One can take for a "constant" 
A^ the value given by Eq. (fT4| , as being the best one 
from all the three, and define a new J-t new as 
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2\Gv\ 2 (l-S R )(l-6 c ) 
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and try to obtain from fit values for 8r and 5n , which can 
be compared with the values computed in [4|. A careful 
analysis of all these numerical values could say that the 
values obtained from fit for 8r and 8c are acceptable, 
and in such a case the problem is closed, or that they are 
not compatible with the theoretical knowledge on T = 1 
analog + states. In the last case one could think to 
small contributions due to scalar and tensor terms in the 



weak interaction model, or that g\/(0) has a small nucleus 
dependence. 

We also did a fit by using directly the values for ft 
given in TABLE IX from [3], such that finally we ob- 
tained four slightly different matrices for the moduli Vij , 
which were used to get a mean value matrix and its cor- 
responding uncertainty matrix. For that we used the 
natural embedment of unitary matrices into the double 
stochastic set, see [l|, or [6l|, for details, and the obtained 
numerical values for moduli Vij and their uncertainties 
are those given in Table. 

Our result on the parameter qaIqv wa s obtained by us- 
ing practically all the measurements where it is involved. 
The values spreading is between 1.25949 corresponding to 
Aq value from paper [53, and 1.27798 obtained by using 
neutron lifetime from Jl2l|. such that the "unexpected" 
Serebrov et al. result, [5J], enters naturally in the game. 

A second conclusion of this approach is that it allows 
a "fine structure analysis" of all experiments measuring 
one definite physical quantity, such as A^, or A, as above, 
providing to each experimental group one measure of how 
far from the ideal situation their measured values are 
standing. 

From fit we have obtained also an experimental value 
on 8 C , which represents the combined radiative and short- 
range physics corrections, see [bH ]. 

The obtained central values for the decay constants 
and form factors are those normally expected, all the 
numbers are given in MeV, and the errors are at la 
level. The results provided for /„. and fx are slightly 
different from that given in (l5j . in particular the lower 
value for /#■. The most critical situation is that for .F(l) 
and G(l), whose central values are almost the same by 
taking into account the huge errors. We remind that 
we used the published results, making no scaling as it 
is usually done, see, e.g., [62( |. The simplest situation is 
that of G(l) where the older data, [4l| and [5l|, give a 
lower value, G(l) = 808.7 ± 77.6, while the new ones, 
[H and [H|, provide G(l) = 1,087.4 ±6.6, which ex- 
plain the huge error. The more complicated case is that 
of .F(l) because of the experimental difficulty to mea- 
sure the product F{l)\U c b\. For example, by selecting the 
highest, and respectively, the lowest values for ^(l), ob- 
tained from data coming from [47| , and respectively [4l[ , 
their ratio is 1.37, which is too big in our opinion. As a 
matter of fact the highest value equals 1, 133.5, which is 
closer to the second value for G(l) given above. Another 
source of error are the assumptions made on form factor 
parameters used to analyse the data, and on the con- 
straints imposed on the shape of these form factors. Of 
course one can exclude some data, but up to now there 
is no definite, or accepted by all, procedure to do it. 

However these problems can be solved in our approach 
since the important measurable quantities from experi- 
ment are \!F 2 (uJi)U c b\ 2 \, and these can be used directly 
in our approach. This procedure can be used to all the 
semileptonic decays, and it will lead to very precise mod- 
uli values, and, more important, to the measurement of 
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form factors moduli in the physical region. It also allows 
a direct test of the lepton universality by putting in the 
fit the measured data obtained separately for electron, 
muon and r, etc. 

More details about our fit, as well as on other physical 
parameters that can be computed by using the above 
results, will be given elsewhere. 



As one final conclusion we can say that by taking prop- 
erly into account the unitarity constraints we found one 
tool which allows the determination of CKM matrix el- 
emnts, and of various decay constants and form factors 
directly from experimental data. 
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